Ultrahigh spatiotemporal Resolution Beam Signal Reconstruction with bunch phase compensation 
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Various electromagnetic signals are excited by the beam in the acceleration and beam-diagnostic elements 
of a particle accelerator. It is important to obtain time-domain waveforms of these signals with high temporal 
resolution for research, such as the study of beam-cavity interactions and bunch-by-bunch parameter measure- 
ments. Therefore, a signal reconstruction algorithm with ultrahigh spatiotemporal resolution and bunch phase 
compensation based on equivalent sampling is proposed in this paper. Compared with traditional equivalent 
sampling, the use of phase compensation and setting the bunch signal zero-crossing point as the time reference 
can construct a more accurate reconstructed signal. The basic principles of the method, simulation, and experi- 
mental comparison are also introduced. Based on the beam test platform of the Shanghai Synchrotron Radiation 
Facility (SSRF) and the method of experimental verification, the factors that affect the reconstructed signal qual- 
ity are analyzed and discussed, including the depth of the sampled data, quantization noise of analog-to-digital 
converter (ADC), beam transverse oscillation, and longitudinal oscillation. The results of the beam experiments 
show that under the user operation conditions of the Shanghai Synchrotron Radiation Facility (SSRF), a beam 
excitation signal with an amplitude uncertainty of 2% can be reconstructed. 
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I. INTRODUCTION 


Many electromagnetic signals are excited by electron 
beams in beam pipes, accelerators, and the beam diagnostic 
elements of particle accelerators [1, 2]. These signals usually 
consist of a mixture of various oscillating modes and change 
rapidly with time [3, 4]. Therefore, obtaining the time do- 
main waveforms of these signals with high temporal resolu- 
tion is important for the study of related problems, such as 
the beam-cavity interaction and the measurement of bunch- 
by-bunch beam parameters. 

These types of signals can be calculated using analyti- 
cal or numerical simulation methods [5-8]. However, it is 
difficult to complete the calculation when the structure of 
the beam vacuum chamber is complex, and it is difficult to 
judge whether the calculation results are sufficiently accurate. 
Therefore, it is necessary to explore experimental methods 
for measuring such signals with high temporal resolution. In 
addition, high-precision bunch-by-bunch measurements have 
become a focus of accelerator research worldwide in recent 
years owing to the advancement of electronic technology and 
the requirement for accelerator machine research [9-19]. 

However, the time resolution (sampling rate) and voltage 
resolution (effective number of bits) of the actual sampling 
equipment are constrained by various conditions [20]. In 
beam measurement, the excitation signal of the electrode is 
often very narrow in the time domain. Taking the Shang- 
hai Synchrotron Radiation Facility (SSRF) as an example, 
the pulse signal of a single-bunch excitation has a width of 
less than 500 ps in the time domain, and the wakefield sig- 
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nal has a width of hundreds of ps [21]. In addition, the am- 
plitude of the wakefield signal is relatively small compared 
to the bunch signal, which requires a larger voltage resolu- 
tion to measure it effectively [22, 23]. Therefore, there is a 
high demand for both time resolution and voltage resolution. 
While conducting this research at the beam instrumentation 
(BI) Group of the Shanghai Synchrotron Radiation Facility 
(SSRF), we could not find any sampling equipment that met 
market requirements. 


A set of technologies for 3D bunch-by-bunch position mea- 
surements based on a high-sampling-rate oscilloscope was re- 
cently developed by the Beam Instrumentation (BI) Group of 
the Shanghai Synchrotron Radiation Facility (SSRF) [24, 25]. 
The core algorithm of this technology is the response function 
projection method of the electrode to bunch signal. To obtain 
ultrahigh time resolution and ultrahigh voltage resolution re- 
sponse functions based on the equivalent sampling principle 
and statistical method, a data processing technique was devel- 
oped to reconstruct the high-resolution periodic beam signal 
using the measured data. 


This paper provides a comprehensive introduction to the 
basic principles, simulations, and an experimental compari- 
son of this method. By conducting beam experiments, the 
relationships between the quality of the reconstructed signal 
and the sampling rate of the data acquisition equipment, the 
depth of the sampled data, and the quantization noise of the 
ADC were analyzed. Additionally, the quality of the recon- 
structed signal under different beam conditions was investi- 
gated, and the impact of transverse and longitudinal beam 
position changes on signal reconstruction was discussed. 


II. BASIC PRINCIPLE 
A. Equivalent Sampling 


To better match the measurement of various parameters in 
the beam measurement, the extraction of the zero-crossing 
point of the pulse signal was used as the time zero point for 
the coordinate axis in this study, instead of using the start- 
ing sampling point, as in traditional equivalent sampling [26]. 
Moreover, this zero point is used to characterize the phase 
zero point of the beam bunch. 

We assume that the excitation signal of the beam bunch is 
sufficiently stable and can be considered as a periodic signal 
with a revolution period within a certain time interval. When 
an ideal data acquisition system is used to digitize the bunch 
signal s(t) with 2N + 1 sampling samples, the result of the n 
th sampling can be expressed as 


x(n) = s(To + nT;). (1) 


where n = —N,—(N — 1),...,—-1,0,1,2...N, To is the 
zero-crossing point of the pulse signal, and 7’, is the sampling 
period. 

Because the signal is considered periodic, it can be ex- 
pressed as 


s(t) = s(t + kT). (2) 


where k = 0,1,2,3,..., T is the period of the signal. 

The same signal is sampled several times by adjusting the 
initial sampling time. For example, by adjusting the initial 
sampling time step to At and the sampling time to M, the 
sampling result can be expressed as 


x(n, m) = s(To + mAt + nT;). (3) 


where n = —N,-(N — 1),.. 
1,2,3...M. 

If the results of multiple samplings are equivalent to the 
same signal period, the sampling rate can be increased. The 
sampling rate is known as the equivalent rate. Theoretically, 
if there is an integer multiple relationship between T, and At, 
the equivalent sampling time interval can be reduced to At, 
which corresponds to the equivalent sampling rate increasing 
from 1/T, to 1/At [27, 28]. 

Figure | (a) shows a typical waveform of the equivalent 
sampled using the dashed blue line. It exhibits a signal with a 
period of 2000 ps. To sample the waveform three times using 
an ideal ADC, 12 points were sampled for each period. Based 
on the principle of equivalent sampling, mapping the data col- 
lected from the three cycles onto a single cycle can result in a 
higher relative sampling rate for the sampled signal, as shown 
in Fig. 1 (b). However, the bunch signal was not strictly 
periodic. The waveform represented by the dashed red line 
in Fig. 1 (a) represents the sampling of a bunch signal with 
longitudinal oscillations. If the waveform is directly recon- 
structed based on the traditional equivalent sampling method, 


.»-1,0,1,2..N,m = 


the results are as shown in Fig. | (c). The blue circles rep- 
resent the equivalent sampling points and the dashed red line 
represents the original signal waveform. It can be observed 
from the Fig. | (c) that the sampling points have deviated sig- 
nificantly from the true values. Furthermore, Fig. 1 (d) repre- 
sents the reconstruction result after applying the turn-by-turn 
bunch phase compensation technique. The use of the turn- 
by-turn bunch-phase compensation leads to a reconstructed 
signal that closely matches the source signal. 
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Fig. 1. (a) Diagram of equivalent sampling (b) Traditional equivalent 
sampling. (c) Traditional equivalent sampling with longitudinal os- 
cillations. (d) Equivalent sampling reconstruction with turn-by-turn 
bunch phase compensation technique. 


B. Signal Reconstruction 


For ideal periodic signals and ADCs, the reconstructed sig- 
nals should exhibit smooth curves. However, actual bunch 
signals are not strictly periodic. It exhibited certain longitudi- 
nal and transverse oscillations, and the ADC exhibited quan- 
tization errors and time jitter. To reduce the error caused by 
beam oscillations, the voltage and time of each sampling (data 


for each ring of the beam) were subtracted from the calcu- 
lated transverse and longitudinal oscillations to approximate 
the real bunch waveform. To reduce the random error intro- 
duced by the ADC, the reconstructed signal was sliced by 
time, and the most probable value in each slice was used to 
represent the value, with the time midpoint used to represent 
the time. 

As shown in Fig. 2 (a), for a signal with a period of 2000 
ps, if an ADC with a time interval of 100 ps (sampling rate 
of 10 GHz) is used for digitalization, 7000 samplings are per- 
formed, with 20 sample points each time and a 3 ps earlier 
starting time for each sampling. After splitting 7000 sets of 
sampling data, a set of reconstructed data with an equiva- 
lent sampling rate of 70THz was obtained. The blue wave- 
form represents the reconstructed signal before subtracting 
the oscillation, and the orange waveform represents the recon- 
structed signal after subtracting the oscillation. In this case, 
the beam oscillation was relatively small, and no significant 
difference was observed, as shown in Fig. 2 (a). In the sub- 
sequent analysis, the impact of the beam oscillation on the 
quality of the reconstructed signal was further discussed and 
experimentally verified. 

The reconstructed signal was divided into 100 time slices, 
with the midpoint representing the time of each slice. Fig- 
ure 2 (c), (d), (e), and (f) show the frequency distribution his- 
tograms of slices 21st, 41st, 61st, and 81st slices, respectively. 
The curves in the figures represent normal distribution fitting 
curves. The value corresponding to the maximum frequency 
was considered the value for each slice. Thus, the response 
function is obtained via spline interpolation as shown in Fig. 2 
(b). 

According to the principle of the equivalent sampling and 
signal reconstruction method, the sampled signal should sat- 
isfy the following conditions: it should be periodic (strictly 
repetitive) or quasi-periodic (with negligible variation during 
sampling) and have a high signal-to-noise ratio. 


HI. UNCERTAINTY ANALYSIS 


In the actual sampling process, the beam signal is not 
strictly periodic and contains nonrandom disturbances, such 
as longitudinal oscillations. These processes are determinis- 
tic and measurable and can be considered during signal re- 
construction to refine the beam signal and make it closer to a 
periodic signal. Additionally, for a non-ideal ADC, the sam- 
pled data contain random noise, which can lead to excessive 
errors in the curve fitting or interpolation results. 


A. Quantizing Noise 


During the sampling process, the quantization noise can be 
approximated to follow a uniform distribution. If the quan- 
tization step size was 6, the quantization noise interval was 
[—6/2, +6/2] [29, 30]. Therefore, the quantization noise can 
be expressed as: 
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Fig. 2. (a) Reconstructed signal of equivalent sampling. (b) Recon- 
structed signal waveform after spline interpolation. (c) Frequency 
distribution histograms of the 21st slice. (d)Frequency distribution 
histograms of the 41st slice. (e) Frequency distribution histograms 
of the 61st slice. (f) Frequency distribution histograms of the 81st 
slice. 


Nq(t) ~ U(—8/2, +8/2). (4) 


where nq (t) is the quantization noise of the sampling point 
at time t and U (—8/2, +8/2) represents a uniform distribu- 


tion over the interval [—5/2, +6/2]. 

Because the quantization noise is additive, when only 
quantization noise is considered, the sampled signal can be 
expressed as 


x(t) = s(t) + n,(t). (5) 


where x(t) is the sampled signal, s(t) is the real signal, and 
n,(t) is quantization noise. 
Therefore, at any time tọ, the sampled signal satisfies 


a(to) ~ U(s(to) — 6/2, s(to) + 4/2). (6) 


Therefore, for time to, the sampled signal of this point sat- 
isfies the probability density distribution, whose mean value 
is the most probable value of this point. 

Based on the above analysis, a quantization noise with 
2% of the original signal amplitude as ô was constructed for 
the simulation. The most probable value for each point was 
considered the result of the reconstructed signal. The rela- 
tive uncertainty of the reconstructed signal was obtained via 
point-by-point subtraction of the original and reconstructed 
signals. The results show that when only quantization noise 
is present in the system, the uncertainty of the reconstructed 
signal is randomly distributed over the entire time domain and 
the relative uncertainty of the reconstructed signal is less than 
1 x 10~? for the quantization noise with ô, which is 2% of 
the signal amplitude. 


B. Time Jitter 


Similar to quantization noise, time jitter is a random signal 
at the sampling time and follows a normal distribution with 
expectation zero and variance o+, which is negatively related 
to the clock stability of the ADC [31]. Therefore, for any 
sampling point, the sampling time tp meets: 


to ~ N(to, ot). (7) 


Assuming that the sampling time to and real signal s sat- 
isfy the mapping relation f, the probability distribution of the 
sampled signal at time tọ is shown in Fig. 3 (a). According 
to the probability distribution of the sampled signal, it can 
be inferred that when the real signal satisfies monotonicity 
near the sampled point, the sampled signal follows a quasi- 
normal distribution, with the real value of the point as the 
most probable value, and the mapping relation f affects the 
sparsity of the probability density on both sides of the most 
probable value. When the real signal has an extreme point 
near the sampling point, if there is f(t,) = f(t2) such that 
p(ti1)+p(t2) > p(to), where p represents the probability den- 
sity at the corresponding time, then the most probable value 
is not equal to the value of the real signal at this point. Oth- 
erwise, the most likely value of this point is the value of the 
real signal at this point. 
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Fig. 3. (a) Probability density distribution of the sampled signal with 
time jitter. (b) Probability density distribution around the negative 
peak. (c) Probability density distribution away from the peak. (d) 
Relative uncertainty of the reconstructed signal. 


Based on the above analysis, the clock jitter with a variance 
of | ps ( the period of the original signal was 2000 ps, and the 
time interval was 0.1ps. ) and mean of 0, constructed. The 
sample number for each point was 10000, and the frequency 
distribution at each point was calculated. Figures 3 (b), (c) 
show the results. As shown in Fig. 3 (b) and (c), the simula- 
tion results are similar to those of the uncertainty analysis. At 
points away from the peak, the sample distribution followed 
an approximately normal distribution, and the most probable 
value was approximately the true value. Near the peak point, 
the frequency distribution of the sample no longer followed 
a normal distribution; however, the most probable value ap- 
proximated the value of the peak point. Therefore, the most 
probable value can also be used to characterize the true value 
of a point. 


The most probable value for each point was considered the 
result of the reconstructed signal. The relative uncertainty of 


the signal reconstructed by the proposed method can be ob- 
tained by subtracting the reconstructed signal from the most 
probable distribution of the original signal, as shown in Fig. 3 
(d). In addition, Fig. 3 (d) shows that the effect of time jit- 
ter on the signal reconstruction is mainly concentrated away 
from the peak, and the uncertainty is small near the peak. For 
a beam signal with a period of 2000 ps and a time interval of 
0.1 ps, the relative uncertainty caused by the time jitter is less 
than 1 x 107°. 

The existing acquisition equipment in our laboratory be- 
longs to Keysight’s Infinitum MXR series, with a clock jitter 
of approximately 120 fs, which is far less than the ps-level 
beam longitudinal oscillation during the user operation condi- 
tions of the Shanghai Synchrotron Radiation Facility (SSRF). 
Therefore, in the actual experiments, the error in the time di- 
mension mainly originates from the beam’s longitudinal os- 
cillation, and the influence of the clock jitter on this method 
can be neglected. 


IV. EVALUATION METHOD 


This paper proposes methods for evaluating the quality of 
reconstructed signals under different conditions. One method 
involves calculating the standard deviation between the re- 
constructed and reference signals, which requires a relatively 
clear reference signal that may not exist in practical experi- 
ments. Therefore, the best measurement result obtained un- 
der the optimal experimental conditions was selected as the 
reference signal to evaluate the relative uncertainty. 

Another method involves evaluating the local quality of the 
reconstructed signal by calculating the second-order differ- 
ence in the signal, which indicates the degree of signal un- 
smoothness. The standard deviation of the second-order dif- 
ference of the reconstructed signal was calculated to evaluate 
the overall quality of the reconstructed signal, which is re- 
ferred to as relative noise in this study. 


V. BEAM EXPERIMENT 


The Shanghai Synchrotron Radiation Facility (SSRF) is a 
third-generation synchrotron radiation source with an elec- 
tron energy of 3.5 GeV, featuring high current intensity, high 
brightness, high beam stability, and low emittance. For the 
storage ring of the SSRF, the RF frequency was 499.654 
MHz, and the harmonic number was 720 [32]. Based on the 
above parameters, to analyze the relationship between differ- 
ent acquisition system parameters, different beam conditions, 
and the quality of the reconstructed signal, a series of beam 
experiments are designed accordingly for verification. 


A. The SNR of Sampled Data 


When the SSRF operates under the user operation condi- 
tions, the filling mode of the bunch in the storage ring is 
shown in Fig. 4 (a). There are 720 buckets in total filled with 


500 bunches, each of which is divided into four bunch trains, 
each containing 125 bunches, with 55 buckets spaced. Each 
bunch train is nearly uniformly distributed; however, there 
are bunches with larger or smaller charges at the head and 
tail of each bunch train. The range is fixed during a single 
sampling event using an oscilloscope, resulting in different 
signal-to-noise ratios (SNRs) between bunches with varying 
charges. Therefore, we evaluated the impact of SNR on the 
quality of the response function by comparing the responses 
of bunches with different charges based on single-sampling 


data collected by the oscilloscope. 
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Fig. 4. (a) Filling pattern. (b) Quality of the reconstructed signal 
with different charges. (c) Relative uncertainty of the reconstructed 
signal with different depths of sampled data. 


The response functions were constructed separately for 
bunches with different charges, and the second-order differ- 
ences of the response functions were then calculated. The 
standard deviation of the second-order difference over the en- 
tire period gives the result shown in Fig. 4 (b), where it can 
be seen that there is an approximately inverse relationship be- 
tween the bunch charge and the relative noise for the range of 
0.1pc to Ipc. Therefore, the SNR of the sampled signal and 
the quality of the reconstructed signal were approximately 
inversely proportional within the experimental measurement 
range. 


B. The Depth of Sampled Data 


In this experiment, the effects of different sampling depths 
were evaluated equivalently by truncating the data for differ- 
ent lengths of the same sampling, where the data length was 
measured in the number of beam turns. From a statistical per- 
spective, the larger the sample size, the higher the expected 
statistical results. Therefore, in this experiment, the maxi- 
mum depth of the oscilloscope was approximately 7000 turns 
selected as the benchmark for calculating the relative uncer- 
tainty of different sampling depths. The results are presented 
in Fig. 4 (c). It can be observed that the greater the number of 
sampling groups, the smaller the relative uncertainty and the 
better the signal reconstruction quality. 


C. Beam Transverse Oscillation 


To assess the impact of beam transverse oscillation on the 
quality of the reconstructed signal, the transient process dur- 
ing storage-ring injection was used as a suitable window be- 
cause it exhibits strong transverse oscillation. Therefore, the 
experimental data selected for analysis contained the storage 
ring injection. As shown in Fig. 5 (a), these data are injected 
after approximately 800 turns, resulting in a transverse oscil- 
lation of approximately +1mm. Figure 5 (b) shows the nor- 
malized amplitudes of the signal peak values obtained by a 
single electrode for the same bunch at different turns, which 
were used to characterize the beam transverse oscillation. As 
shown in Fig. 5 (b), there was a relative change of 10% after 
the injection. 

Signals were reconstructed using data from the first 800 
turns during injection. The first 800 turns were not affected by 
the beam injection; therefore, they were used as benchmarks. 
The relative uncertainty was obtained by subtracting the two 
reconstructed signals, as shown in Fig. 5 (c). From these re- 
sults, it can be speculated that the relative uncertainty of the 
reconstructed signal can reach approximately 2% when there 
is a 10% transverse oscillation in the beam signal. Therefore, 
under the condition of transverse oscillation, if the required 
accuracy of the constructed signal is less than 2%, the influ- 
ence of transverse oscillation can be ignored, and the signal 
can be directly constructed; however, when the required ac- 
curacy is better than 2%, the data with small transverse oscil- 
lation or the refined method should be selected to construct 
the reconstructed signal. 


D. Beam Longitudinal Oscillation 


To analyze the influence of longitudinal oscillations on 
the quality of the reconstructed signal, beam signals contain- 
ing longitudinal oscillations were selected for evaluation, as 
shown in Fig. 6 (a). The longitudinal oscillation ranges from 
-6 to 2 ps. 

Based on the actual longitudinal oscillation, the original 
data were refined by subtracting the bunch phase caused by 
the longitudinal oscillation per turn, and the reconstructed 
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Fig. 5. (a) Transverse position of different turns. (b) Amplitude 
of the signal of different turns. (c) Relative uncertainty induced by 
transverse oscillation. 


signal was constructed based on equivalent sampling. The 
difference between the refined results and the results of direct 
signal reconstruction without compensation is the relative un- 
certainty, as shown in Fig. 6 (b). As the reconstructed signal 
was normalized, the result in Fig. 6 (b) can also be viewed as 
the relative uncertainty from the longitudinal oscillation, with 
the relative uncertainty being better than 2 x 1078. 

Furthermore, the second-order difference in the recon- 
structed signal before and after compensation can be calcu- 
lated as relative noise, as shown in Fig. 6 (c) and (d). Com- 
pared with the reconstructed signal, it can be seen that the 
relative noise is larger at the point with a larger slope of the 
reconstructed signal, and this part of the difference can be 
significantly reduced by refining the original data. 

In addition, the storage ring may be in different longitu- 
dinal oscillation states according to the machine conditions, 
and another set of data with different oscillation states was 
selected for the experimental analysis. As shown in Fig. 6 
(e), the longitudinal oscillations range from a few ps to -300 
ps. Similarly, the results in Fig. 6 (f) can be obtained by re- 
constructing the original signals before and after compensa- 
tion and then taking the difference as the relative uncertainty. 
Clearly, when the amplitude of the longitudinal oscillation is 
large, the uncertainty is very large, and the relative uncer- 
tainty reaches 10% under the conditions of this experiment. 

In summary, for the different acquisition systems, the 
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Fig. 6. (a) Variation diagram of longitudinal oscillation. (b) Rela- 
tive uncertainty from longitudinal oscillation. (c) Relative noise of 
the reconstruction signal before compensation. (d) Relative noise of 
the reconstruction signal after compensation. (e) Variation diagram 
of the larger longitudinal oscillation. (f) Relative signal difference 
between amplitude before and after compensation 


higher the number of ADC sampling bits and the larger the 
sampling depth, the better the quality of the reconstructed sig- 
nal. For different beam conditions, longitudinal oscillations 
with large amplitudes must be compensated before signal re- 
construction; otherwise, the authenticity of the reconstructed 


signal cannot be guaranteed. However, for transverse and 
small longitudinal oscillations, compensation can be selected 
according to the system requirements. 


VI. PERFORMANCE EVALUATION 


Based on the beam experiment, 200 sets of beam data 
at different times with a maximum sampling depth (7000 
turns), maximum sampling signal-to-noise ratio, and negli- 
gible transverse and longitudinal oscillations were selected 
for signal reconstruction. The standard deviations of 200 
sets of reconstructed signals were calculated point-by-point 
after normalization. The relative amplitude uncertainty was 
characterized by this standard deviation, as shown in Fig. 7. 
Therefore, for the signal reconstruction method proposed in 
this study, the relative amplitude uncertainty is better than 2% 
when the sampling depth is at its maximum, the sampling 
signal-to-noise ratio is at its maximum, and the lateral and 
longitudinal oscillations are negligible. 
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Fig. 7. Relative uncertainty of reconstructed signal at different times. 


VII. CONCLUSION 


A reconstruction method for periodic beam signals in an 
electron storage ring was developed based on the principle of 
equivalent sampling. Compared with traditional equivalent 
sampling, when constructing the reconstructed signals of each 
bunch, the effect of multiturn bunch instability is reduced by 
subtracting the transverse and longitudinal oscillations. The 
zero-crossing point of the signal pulse is taken as the time 
zero point, which is conducive to the subsequent study of sig- 
nals, such as the wakefield of the beam based on the bunch. 

Beam excitation signals with a relative amplitude uncer- 
tainty of 2% could fully satisfy the requirements of the re- 
sponse function of the bunch-by-bunch 3D measurement sys- 
tem. The reconstructed signal bandwidth in this paper was 
4.2 GHz, and this method has no bandwidth limitations. If 
the data acquisition equipment and signal transmission net- 
works are upgraded, they can be used to process signals with 
higher bandwidths. Therefore, this method can be extended 
to the sampling analysis of wakefield signals with bandwidths 
of tens of GHz bandwidth. 
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